Non-Equilibrium Duality 
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The dual 'worldline', or 'string', description of adiabatic Ginzburg-Landau field theory near a 
phase transition in terms of quasi-Brownian strings and loops is well understood. In reality, the 
implementation of a transition is intrinsically non-equilibrium. We sketch how time-dependent 
Ginzburg-Landau theory leads to a modified dual string picture in which a causal bound on the 
growth of unstable string prevents the uncontrolled proliferation of string (the Shockley-Hagedorn 
transition) suggested by the adiabatic approximation. 
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From the early days of lattice models jjj it has been 
appreciated that time-independent Ginzburg-Landau 
(TIGL) theory has a dual 'worldline' representation in 
terms of near-Brownian strings and loops. More details 
are given in Refs. [0,||. 

The main ingredient in all these calculations is that 
they are performed for systems in equilibrium, in which 
the temperature T is fixed at values close to the critical 
temperature T c . For the case that interests us, that of 
continuous transitions, the instabilities that characterise 
the transitions are, in this adiabatic dual picture, a con- 
sequence of the uncontrolled proliferation of string. 

In practice, transitions occur in a finite, often short, 
time in which the temperature T crosses T c from its ini- 
tial value. The adiabatic approximation is only valid 
away from the transition, when the field can keep in step 
with the changing environment. Close to the transition 
the adiabatic regime is replaced by an impulse regime, 
in which the field falls out of step with what its equi- 
librium behaviour would be, and freezes in ||||. Just 
as the correlation length, in reality, is unable to diverge 
because there is not enough time, so we would expect a 
realistic dual worldline description in terms of strings to 
show a limit to their proliferation as we pass through a 
transition. 

The question that we shall begin to address here is 
how the intrinsically non- adiabatic transition has a dual 
representation in terms of strings and loops. 

However, we begin with a brief recapitulation of equi- 
librium dual theory. The simplest condensed matter 
TIGL theory is that of a single complex field (f>, with 
Ginzburg-Landau free energy 



F(T) = J d 3 x (^L\V^ + a{TM< 



in which the chemical potential a(T) = a^e{T), where 
e(T) = (T/T c — 1), vanishes at T c . Such an energy pro- 
vides a reasonable description of superfluid 4 He, a sim- 
plified model for 3 He and a good description of the scalar 
sector of low-T c superconductors. 

It is convenient to work in spatial units of £o = 



h 2 /2mao when, on rescaling the field, self-coupling and 
temperature, 

F(T) = / d 3 x (|V</>| 2 + e(T)|0| 2 +^| 4 



At temperature To > T c the free-field correlation function 
that follows from 



is 



F (T )= / d 3 x |V0| 2 + e(To)H 



(r=| 



Wx)f(0)> = G (r) = J fke^P(k), 
x|) in which the power spectrum 



P(k) 



1 



k 2 + e(T ) 



dre- Tk \~ r ^ 



has the usual representation in terms of the Schwinger 
proper-time r. In turn, this gives 



Go(r) 



dr[ 1 
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Go(0), necessary for loops, diverges from the UV singu- 
larities at t — 0, but these can be regulated with a fixed 
cutoff. We choose a cutoff of order unity in units of t;o, 
which will be implicit throughout. 

The dual picture is obtained by observing that 



(1/47TT) 



3/2„-r74r 



is the probability distribution for a 



Brownian 'worldline' or, more usefully, a 'polymer', path 
x(t), parametrised by r, having its endpoints separated 
by a distance r. Specifically, it is expressible as the sum 
over such paths 



3/2 



f-V • 

\4ttt J 



■r 2 /4r 



x(t)=x 



£>x exp 



(0)=o 



If we think of the path as having step length unity (in 
units of £o) then r is proportional to the length of the 



1 



path H and we shall use r and path length synony- 
mously. 

Then Go(r) is the sum over string paths of all lengths, 

poo />x(t)= X 

G (r) = dT p xe -S e5 [x i r,e(T )] ! 

JO Jx(0)=O 

where S^x; r, e(To)] is the equilibrium Euclidean action 
5 e9 [x; r, e(T )] = jf dr' ( \{^f + e(T )) . 
The free-field partition function 
Z= [ V<j)V(t)* e- p ° = exp[-<rln(-V 2 +e(T ))] 



is just that of a gas of free orientable polymer or string 
loops, whose lengths are parametrised by r, with one-loop 
partition function 



InZ = 



dr ( 1 
t \4ttt 



3/2 



-Te(To) 



(2) 



e(To) is understood as the energy/unit length or tension 
of the string. The factor oc t~ 5 / 2 in (|^) is the length 
distribution for Brownian strings of length r. As we 
drive To — > T c + it is easier to create string. If, formally, 
we take To below T c then the onset of 'negative 1 ten- 
sion makes long loops overwhelmingly favoured and the 
condensation of loops that follows is understood by con- 
densed matter physicists as a Shockley-Feynman transi- 
tion, and by quantum field theorists as a Hagedorn tran- 
sition. 

The average loop length is 



r ) = - 



JO \ 4-7TT 
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e(To) 



d In Z 

de(T ) ~ / \ 3 / 2 

JO T \ iTTT 1 e 



(3) 



Yet again, forcing T below T c makes (r) diverge from 
the IR divergence of the integrals at large r. A formal 
IR cutoff at t = T rnax leads_to (r) = 0(r max ). 

The incorporation of the f3\(f>\ field self-interaction into 
equilibrium dual string theory is implemented by the in- 
troduction of a repulsive steric string interaction at the 
points where paths or strings cross ||. Although this 
reduces the weight of string configurations, qualitatively 
the transition still occurs in this adiabatic picture be- 
cause of the proliferation of strings. 

In practice, a change of phase is enforced by chang- 
ing the sign of e(t), either by reducing T, or by varying 
the critical temperature T c . Experiments with supercon- 
ductors and (the neutron bombardment of) 3 He do the 
former, while pressure quenches of A He do the latter. If 
we write 



e(t) = e(T(t)) 



T(t) 
Tc(t) 



- 1, 



the adiabatic approximation for the 'free' correlation 
function Go(r,t) at time t is 

poo px.(r)—x 

Gg d (r,t)=/ dr Dxe^^'O', (4) 

JO Jx(0)=O 

in which we just make a straightforward substitution of 
the equilibrium c(Tq) with e(i) in 5 eg [x, r, e(T )]. This is 
treating the equilibrium pictures as a series of snapshots 
that can be run together as a continuous film. In this 
film string lengths increase uncontrollably as we cross 
the transition. In particular, Gg d (r, t) of (|l|) is simply 
calculable as the Yukawa correlator 



Gf{r) 



where, on rescaling, £ ac t(t) — fo/\AW diverges as 
T(t) — » T c . This cannot be the case, since causality alone 
prevents the correlation length diverging in a finite time. 
In terms of the dual picture this implies that the pro- 
duction of an infinity of string in a finite time is equally 
prohibited. The question is, how? 

To answer this, we adopt the time-dependent Landau- 
Ginzburg (TDLG) equation for F, 



1 d(t> 5F 



V, 



(5) 



r dt 5<p 

where r\ is Gaussian thermal noise, satisfying 

fo(x, t)rf (y, 0> = 2T(t)r<5(x - y)S(t - t'). (6) 

Let us continue to consider the free-field case, F = /Jo- 
in (J^j the natural unit of time is To = 1/aoT and, in units 
of tq and £oj Eq.(^j) becomes 



(x, t) = -[-V 2 + e(i)Mx, t) + t?(x, t). 



(7) 



where r\ is the renormalised noise. It is straightforward 
to show that the equal-time correlation function is now 

(^(x,t)^(0,t)) = G (r,t) = J fke^P(k,t). 

in which the power spectrum P(k, t) has a representation 
in terms of the Schwinger proper-time t as 

P{k, t)= dr f(t - t/2) e~ Tk2 e- Jo dT < l - T ^ , 
Jo 

where T is the renormalised temperature. In turn, this 
gives 

GoM) = 

i / 1 \ 3/2 rT 

drf(t-r/2) ( —) e- r2 ^ e ~f ^ 



\4tVT 



(8) 
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This differs significantly from (Q), which would replace 
e(t - t' 12) in (|)with e(t). Further, because of f(t - 
t/2), the paths are no longer Brownian in their length 
distribution. 

However, Go(r, i) can still be expressed in terms of 
paths, as 

/>oo />x(r)— X 

G (r,t) = d T f{t-T/2) Dxe- s[x ' T ' f| . 

JO ix(0)=O 

S^x, r, t] is not the equilibrium action S eq {x., t, e(t)], but 



(i<^» 2+ '< f -^ 2 >)- 



5[x, r, t] 



Unlike the case for the adiabatic approximation (Q), the 
local tension e(t — t/2) is dependent on 6o£/i the time t 
at which we examine the strings and on the position r 
along them. As a result the integrated tension JJ" dr'e{t— 
t'/2) does not vanish at the transition time. Since it 
is the uniform vanishing of tension which triggers the 
avalanche of string production we see already that it will 
not happen in this case. 

As a simple, if unrealistic, demonstration that transi- 
tions implemented in a finite time do not impose singu- 
lar dual behaviour we consider the case of an instanta- 
neous quench at time t = from a temperature To above 
the critical temperature T c to absolute zero. That is, 
e(t) = e(To)0(-t). 

Simple calculation shows that Go (r, t) takes the form 
of (Q) for t < 0, whereas for t > 

G (r,t) = e 2tTo ^ / dr / p xe -S e<! [x,T,e(To)] (g) 

3/2 



/■oo / 1 \ 



(To) 



After the transition we have a representation in terms of 
positive tension paths at the initial temperature To- The 
unstable (negative tension) paths, for which r < 2t, are 
totally excluded. Instead, the instabilities are encoded in 
the non-singular exponential prefactor. 

This lack of IR singular behaviour is made even more 
explicit in a saddle-point approximation for Go(r, t). For 



> 4e(T )i we find 



GoM)~e 2tT ^ — 



-r/t 



Airr 



where £ = £(T)) remains frozen in at its initial equilib- 
rium value. Unlike £ a d(t) there is no divergence of £ as 
we cross T = T c in this abrupt way. 

For the equilibrium theory, (r) of (^) can also be ex- 
pressed as 



(r) 



Gq(0) 
2G '(0)' 



where the prime denotes differentiation with respect to 
r. If we adopt the same definition out of equilibrium the 
benign effect of the quench to Tf = is again apparent 
in that (r) t agrees with (||) for t < 0, but is 



(r)t 



3/2 



r°° dr_ I 1 

Jit T \ 47T 



3/2 



(10) 



for t > 0. The exponential prefactors have cancelled to 
reproduce the equilibrium result @), again at the initial 
temperature T — To, but for the absence of unstable 
loops with length less than 2t, to which there is now no 
reference. Because of the exponential damping of long 
string the dominant string length is r = 2t. 

There is no question of an IR divergence of loop length 
as naively suggested by the equilibrium theory. We un- 
derstand the stability of loops with length r > 2t in ( |To| ) 
as a causal bound. In our units, negative tension can 
only propagate at speed c — 1, which happens to be the 
cold speed of sound in the <j> field. 

This is reinforced by an extension to non-zero final 
temperature Tf 

e(t) = e(T ) > if t < 0, 
= e(T f ) < if t > 0. 

For t > we now find 

/>oo />x(t)— X 

Go(r,t) = e 2t ( T °- T /)/ T = / dr Pxe -s e5 [x,r, e (T o) ] 

rri fit /»x(r)— X 

To Jo Jx(0)=O 

We see explicitly how the limited length r < 2t of string 
with negative tension e{Tf) < prevents such string from 
giving a divergent contribution. Small unstable loops are 
no longer precluded, but their contribution is finite. The 
dominant length remains r = 2t. 

Although an instantaneous quench is impossible, more 
general quenches show similar qualitative behaviour. 
Specifically, suppose that e(t) decreases monotonically, 
with a single zero e(0) = at time t = 0. As before, only 
strings of limited length t < 2t, for which e(t — t/2) < 0, 
have negative tension. However, for r > 2t, for which 
e(t — t/2) > 0, strings have segments with negative and 
positive tension. There is a causal bound c = 1 on the 
speed at which instability can propagate along a string. 
The string with negative tension, with a contribution 
that is independent of t, gives a prefactor growing at 
least exponentially. The effect is to leave only stable 
string of length t > 2t in the t integral for r > 2t. This 
boundedness on unstable string prevents the unlimited 
production of string suggested by the adiabatic approxi- 
mation. 
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The exponential growth of Go(r,t) with time in (10), 
and more generally, can only be accommodated for a very 
short period, since (\(f>\ 2 ) = Go(0, t) must be constrained 
by the value of the order parameter (\<p\) after the tran- 
sition, equal to y /3 _1 /2 in the absence of corrections. A 
rough guide to the maximum time t for which the free- 
field approximation is valid is that Go(0,i) = /3 _1 /2. As 
t approaches t, and thereafter, the reaction of the field 
with itself will cut off the exponential growth. 

It is difficult to see how the dual string picture will sur- 
vive at later times without further approximation. One 
indication is through a mean-field (or large-N) approach 
in which the linear nature of the theory is maintained 
self-consistently || . In this approach (Q) is replaced by 



(x, t) = -[-V 2 + e e// (i)]0(x, t) + f/(x, t). (11) 



with 



eff(t) 



e(t) 
e(t) 



where G(0, t) is determined self-consistently from (11). 
The coefficient p depends on whether we adopt a mean- 
field (Hartree) approximation or a large-N limit for N = 
2. 

The assumed single zero of e(i) at t = will lead to 
a zero of e e //(i) at t w 0, and the previous analysis ap- 
plies; with only segments no longer than 2t with negative 
tension, there is no singular behaviour as we cross the 
transition. 

To have a quantitative estimate of the effect of back- 
reaction we make a Gaussian approximation for G(0, t) 
by expanding about the zero of e e //(i). Assuming that 
the quench is not too rapid, we find that, for t < i, the 
average loop length is 



(r)t 



roc , f(t-r/2) -(r-2t) 2 |e'(0)|/4 
JO (A-kt) 3 / 2 C 

Jo 



(4;tt) 3 / 
T(t-r/2) 
(4ttt) 3 / 2 C 



(r-2t) 2 | e '(0)|/4' 



As before, the prefactors encoding the unstable strings 
cancel approximately, and we need no information about 
the self-consistent mass, but for the fact that e e //(i) — > 
at large times to stop G(0,i) from growing. We see the 
peak at r = 2t in the length distribution moving clear 
of the UV endpoint behaviour. This continues for later 
times and (r) t = 0{t) once the temperature has become 
low enough that T(t) has suppressed the UV singular 
behaviour. Details are given in our earlier work (7|^]. 

We stress that this discussion has been restricted to 
the 'first quantised' dual representation of the Ginzburg- 
Landau theory. However, it is a familiar result from a 
different viewpoint. For a linear system like (|i"i"l ) it can 
be shown |lQ,Ol that the mean loop length (r)t satisfies 



where n(t) is the density of line zeroes of the complex field 
cj). The relevance of this is that, at later times, the global 
vortices of this U(l) scalar theory can be identified by the 
line zeroes of their cores Jl2). The linear growth of dual 
loop lengths with time corresponds to a t 1 ' 2 behaviour for 
cj)-field line-zero separation and, when vortices are well- 
defined, vortex separation ]l2| |. This scaling behaviour 
can be justified [[13| for the decay of vortices of 4 He, and 
is used to determine initial vortex densities at the 4 He 
transition p4|] . 

This work is the result of a network supported by the 
European Science Foundation. I thank Adriaan Schakel 
for helpful comments. 
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